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We demonstrate the appearance of zero-energy bound states satisfying the Z 4 parafermionic algebra in inter¬ 
acting quantum wires with Rashba spin-orbit coupling and proximity-induced superconductivity. The fourfold 
degeneracy of these states is protected by time-reversal symmetry and fermion-parity conservation, and leads 
to an 8 ;r periodicity of the Josephson current due to the tunneling of fractionalized excitations with charge e/2. 
Even in the presence of perturbations, we propose that this periodicity will remain visible in driven, current- 
biased Shapiro step measurements on current state-of-the-art Rashba wires. 

PACS numbers: 71.I0.Pm, 74.45.+c, 05.30.Pr 


I. INTRODUCTION. 

One of the most active fields in physics in recent years has been the search for Majorana fermions GM- Despite being 
theoretically predicted more than 70 years ago 0, the existence of Majorana fermions as elementary particles is still unclear. 
However, it was recently recognized that they can occur as quasiparticles in various solid-state systems IMIl- For instance, 
certain one-dimensional (ID) wires are believed to host localized Majorana bound states (MBS) due to the combined effects of 
Rashba spin-orbit coupling (RSOC), proximity-coupling to a superconductor, and an externally applied magnetic field lfT3lfT4l . 
These MBS can be understood using simple single-particle quantum mechanics. Exotic and closely related bound states, called 
parafermions, have been predicted to descend from strongly correlated phases such as the edges of fractional quantum Hall 
systems MM- 

The practical interest in these unusual quasiparticles stems primarily from their suitability for topologically-protected quan¬ 
tum computation IT]. Exploiting their non-Abelian exchange statistics, many logic gates can be implemented by braiding, an 
operation which is robust to local perturbations. Indeed, the prospect of intrinsically decoherence-free qubit operations is a 
strong driving force for this research field, and has the potential to revolutionize quantum computation. 

If experimentally confirmed, MBS and parafermions will be the first discovered particles with non-Abelian exchange statis¬ 
tics. Edge states of two-dimensional topological insulators, ferromagnetic chains on superconductors, as well as semiconductor 
nanowires with RSOC have recently shown experimental signatures that are consistent with the appearance of MBS lfT9l422l . 
In particular, recent work Il23lj25l suggests that it is possible to observe An periodic components of the Josephson current aris¬ 
ing from Majorana fermions, giving a window on their twofold degeneracy. Measured current-voltage curves for proximitized 
edge states of HgTe and nanowires with RSOC have shown the disappearance of odd numbered “Shapiro steps” consistent with 
theoretical predictions In contrast, experimental results indicating the presence of parafermions do not yet exist, largely 
because the proposed experimental setups are difficult to realize. 

In this paper, we aim to tackle this lack of experimental realization of parafermions by describing a rather simple arrangement 
in which a single interacting ID Rashba wire exhibits fourfold degenerate zero-energy states obeying a Z4 parafermionic algebra 
(see Refs. 12614^ and references therein). Our proposed system exploits a generic interaction process which has so far not 
been examined in quantum wires, called spin-umklapp scattering, whereby two spin-up fermions are scattered into two spin- 
down fermions. Such an interaction arises from virtual transitions between subbands, and for a suitable choice of the chemical 
potential opens a partial gap in the energy spectrum of the wire. Adding proximity-induced superconductivity then leads to the 
emergence of parafermion bound states at the ends of the wire, whose fourfold degeneracy is protected by Kramers’ theorem as 
long as time-reversal symmetry (TRS) is not broken. As a definitive experimental signature, we point out that tunneling of the 
fractionalized quasiparticles with charge e/2 will result in an Stt periodic component of the Josephson current, which could be 
seen in Shapiro step experiments. We note that the required interaction strengths have already been reached in quantum wires 
El, and that the opening of the spin-umklapp gap may already have been seen in recent experiments (e.g. M)- Moreover, 
the predicted fourfold degeneracy and the associated Josephson effect are robust to weak disorder, of the magnitude of the gap 
which is opened by spin-umklapp scattering. 
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FIG. 1: A process in which virtual transitions (g) between the lowest and the first excited subbands of the quantum wire allow for a spin-non- 
conserving umklapp process to be generated from a spin conserving interaction V. 


In section 2, we introduce the model for an interacting Rashba wire and derive the spin-umklapp scattering term. We then 
bosonize this model in section 3, and present a renormalization group (RG) analysis for the flow of the system parameters. We 
find there exist regimes in which fourfold degenerate edge states can be found. In section 4, we describe the implications of 
the degenerate ground state for Josephson effect measurements, and suggest a definitive Shapiro step measurement to identify 
the 8;7r contribution to the Josephson current arising from parafermions, even when this contribution is sub-dominant. Finally, in 
section 5, we discuss the relevance of our work for experimental investigations. 


II. MODEL. 


We consider a long quasi-lD nanowire along the x direction which is harmonically confined in the y and z directions. The 
interplay of the intrinsic spin-orbit coupling of a material and the breaking of inversion symmetry of a particular geometric 
arrangement, due to either the presence of a substrate, or to the application of an out of plane electric field, gives rise to RSOC 
ED- In the latter case, the strength of the electric field may be used to tune the magnitude of the RSOC, denoted as or. The 
dynamics in the z-direction is not affected by the Rashba coupling and can be safely ignored. Thus, we can model a finite-width 
wire by using the following 2D Hamiltonian including RSOC with strength a* Il32l434ll 


H = 


•Fy 


2m 


2 2 
-maj y ■ 


■ auicr^py - a-yPx), 


( 1 ) 


where p^ and py are the momentum components in the x- and y-directions, and (Tx^y are Pauli matrices. The transversal confine¬ 
ment is modelled as a harmonic potential with frequency u. The system has translational invariance along the .r-direction but is 
strongly confined along the y-direction, which leads to the appearance of higher excited bands separated from the lowest band 
by a spacing determined by the inverse width of the wire. Introducing raising and lowering operators, a ’ and a, one finds that 
H - Hq + Hi, where 


Ho = + 5 ) + - aRO-yPx, ( 2 ) 

Hi = iga-xia^ - a), (3) 


and g - Ur yfmcjfl. Since eigenstates of Hq have a quantized spin in cry direction, the form of Hi ensures that transitions 
between neighboring subbands, which are caused by a and a *, always involve a spin flip. We will account for the possibility of 
virtual transitions between the lowest subband and the first excited subband with energy u by integrating out the coupling Hi 
between the subbands. We do this by means of a Schrieffer-Wolff transformation, H' - He^ where [5,//o] = Hi. Up to 
second order in ma^joj one finds 

, 1 p\ 

H ^ Ho - -[5,Hi] = -- aRO-yPx, (4) 

2 2m 

which tells us that it would be consistent to simply ignore the dynamics in y-direction in Eq. 0 m. Deviation of the spectra 
from the parabolic form occurs only at order (ma^ju)^ Il32l[^ . and is not relevant for our current discussion. Note, however, 
that the associated change in Fermi velocities can have interesting effects for weak interactions llTSll . 

In the interacting case, the interaction term also gets modified by the Schrieffer-Wolff transformation. A generic density- 
density interaction potential y(ri - r 2 ) where r, = (x,,y,), clearly conserves spin. However, since virtual transitions to higher 
bands come with a spin-flip and the interaction potential y(r) mixes states in different subbands, spin-umklapp scattering is 
generated as shown in Fig. 

After the Schrieffer-Wolff transformation, the interaction potential for the lowest subband takes the form V' - e Ve^ - 
Tp + Fsf + Vsx, with a density-density interaction Vp, a spin-flip term Fjf which changes the total spin of the interacting particles, 
and a spin-exchange term Fjx, which allows particles in the lowest band to exchange their spins, but conserves the total spin. 
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FIG. 2: (color online) This figure shows the position of the chemical potential and the four low-energy linearized modes used to bosonize the 
system. Arrows indicate the RG-relevant interaction processes between these modes: umklapp scattering between the two modes near fc = 0 
(dotted lines) and proximity-induced j-wave superconductivity, which pairs modes of opposite physical spin (double line). 


Placing the chemical potential at the Dirac point, as shown in Fig.|^ gives four low-energy modes at the momenta k -0, ±kf, 
where kf - 2maR. Correspondingly, for small energies, we can split the field operators up into four modes. 


^ Mix) + e-^’^'^Mlix)- 

Next, we express ^asix) for a e [L, R} and s e {"f, i) in terms of its Fourier components, 

^as(x) — ~ ^ C 4^as,k'> 

yL t. 


(5) 

(6) 


(7) 


where L is the length of the wire. In terms of these operators, the interaction Hamiltonians after projection to the lowest subband 
can be written as follows. The density-density interaction term reads 



dxpas(x)PaAx), 


( 8 ) 


where the fermionic densities are defined as usual as Pas(x) — i)/ts(x)i)/as(x) and V(q) is the Fourier transform of the interaction 
potential y(r) projected to the lowest subband. 

Due to momentum conservation, the spin-flip terms only mix terms near k -0, 


Vsf = Vsf J dx + h.c.], 


(9) 


where we retained only the leading local part Vsf = -2t/(0) oc matlcj^, and U(q) is given in Eq. (A8 i. The full calculational 


details can be found in Appendix A. It is the Hamiltonian Vsf which allows for spin-umklapp scattering. 

Most of the spin-exchange terms in Vsx can be expressed as density-density interactions, leading merely to changes in the 
coefficients of the terms in Eq. We separate out the single non-density-density term 


Vs = vs J + h.c.j, 


( 10 ) 


which corresponds to an interaction between inner and outer bands with strength vs = 2kj^U{kf). In the limit in which this term 
dominates, it results in a spin-density wave state atq- 2kf, and has been discussed in detail in Ref. Il34ll . 

Einally, we allow for the possibility of proximity-induced coupling to an s-wave superconductor, which pairs spin-up and 
spin-down electrons, and so has an effect at all the Eermi points in Eig|^ This pairing contribution to the Hamiltonian is 


Vsc^vsc J dx(^l^^l + ^l>/rl^+h.c.), 


( 11 ) 


where vsc is determined by the strength of the proximity-coupling to the superconductor. We now analyse the competition 
between these three possible interaction channels, parametrised by Vsf, vs and vsc- 
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III. BOSONIZATION & RENORMALIZATION GROUP (RG) ANALYSIS. 


To further analyze the interacting system, we write the Hamiltonian in terms of bosonic operators 4>± and 6± by defining 




^Lj 

y/lna 






( 12 ) 


where rjaa- are Klein factors and a denotes the short-distance cutoff. Here, 4>+(x) and 0+(.x) are canonically conjugate bosonic 
operators for degrees of freedom near k = ±kf, whereas (p-ix) and 6-(x) describe modes near k - 0. In these variables, the 
Hamiltonian consists of two Luttinger Hamiltonians for the -H and - species, with approximately equal Luttinger parameters 
K±, and interaction terms reflecting Eqs. (|^ and ( [T0| ). Moreover, one obtains derivative terms which couple the two species 
and Kgd^0+d_,0^. 

Following Ref. Il34l . we can diagonalize the quadratic parts of the Hamiltonian by going to the charge-spin basis (pp^a- — 
{([>+ + 4>S)I V2 and = (0+ + 0_)/ V2 so that 


Hq 



(dxpaf 

Ka 


+ Ka(d,0af 


(13) 


where Vp o- are the respective sound velocities of the modes. For repulsive interactions, we have Kp < I and Ku- < 1 El- In 
addition to Hq, we obtain two competing interaction terms. 


^ J dxcos[2y/20^], (14) 

Vsf ^ J dx cos[2y/2((f)p-(pa-)]- (15) 

Proximity-induced i-wave superconductivity gives a contribution which reads in bosonized form, 

= / dx{cos[^/2(0p + 0^)] + {0^^-0^}'j. (16) 

For weak interactions, all parameters of the model can be determined precisely in the bosonization procedure (see Appendix 
A). However, for strong interactions it is more convenient to regard the parameters Vp_o- and Kpa- as well as the three coupling 
strengths gs, gu 8sc effective parameters, which may flow independently under renormalization as we change the cut-off 
a. 

We calculate the flow of the various coupling constants using real-space RG calculation based on operator product expansions 
ES). We And the following first-order RG equations for the coupling constants of the cosine terms (see Appendix B for details 
of the RG proceedure). 


dgs 

d( 

00 

1 

(N 

II 

(17) 

dgu 

dt 

1 

1 

1 

cT 

II 

(18) 

dgsc 

di 


(19) 


implying that the spin-density wave term is always irrelevant for repulsive interactions (Ka- < 1) ED- The spin-umklapp term, 
by contrast, can become relevant for strong interactions where Kp H- Kg- < 1. Finally, the superconducting term is relevant for 

Kp^+Ka^<4. 

We would like to point out that for = 0, the system becomes SU(2) invariant. In that case, the spin-umklapp term vanishes 
because spin is conserved and one finds the well-known Kosterlitz-Thouless RG flow which brings Kg- 1 as gs ^ 0. In 
contrast, for qr + 0, Kg is not constrained and strong repulsive interactions lead to Kg <si 1. 

We start the RG flow from an initial value a — ao and flow towards a ~ L, the length of the wire. Generically, the RG flow 
will stop at a finite value Qoa < Las soon as one of the dimensionless coupling constants gsc.u approaches one. The bare value of 
gsc(ao) is determined by the strength of the proximity coupling to the superconductor, which can be experimentally optimized. 
The bare value gu(ao) depends on the separation between the lowest subbands, and so depends on the transverse confinement 
(i.e. the physical width) of the wire. 

To generate zero-energy bound states, spin-umklapp scattering must gap out the modes near k = 0, whereas proximity-induced 
superconductivity should open a gap for the modes at k = +kp, see Fig. Superconductivity affects all modes, so this is only 
possible if at the end of the RG flow |gu(aoo)| > lgsc(‘5ioo)| > 0. Strong electron-electron interactions result in Kp < 1 /2 and 
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Ka- < I /2, which a priori makes the spin-umklapp term relevant and the superconducting term irrelevant. However, since the 
RG flow is cut off at a finite length scale, one will generally find a nonzero |gsc(floo)| > 0 at the end of the RG flow, meaning that 
a superconducting gap will still open. 

To demonstrate the existence of degenerate, localized zero-energy bound states, we use the unfolding transformation described 
in Refs. 1221ES). This transformation can be used to map our system with length L and open boundary conditions to a system of 
length 2L and periodic boundary conditions. Explicitly, we construct the unfolded chiral fields 


i: i VR+ix) for 0 < i < L 

^ \ >fL-i2L -X) for L < x < 2L 

e i VL+ix) for 0 < X < L 

’ \ (Pr-(2L - x) for L < X < 2L 


( 20 ) 

( 21 ) 


where ipav - ctcpy - 6y with a - R,L and v = +. In order that the fermionic fields obey the vanishing boundary conditions at 
X = 0 and x = L, we find that the bosonic fields must obey 

tpL+(0) = (fiR-iO), tpL+iL) = tPR-(L) 

tpL-iO) = 9?^+(0), tpi-iL) = </?fi+(T) (22) 


Note that in this transformation, the degrees of freedom 0+) are mapped on the range x e [0, L], whereas the ((^_, 0_) fields 
are mapped on the range x e [L, 2L]. Since the original chiral fields satisfy {ipav{x),ipa’v'ix')} - inadaa'^w’Sgnix - x'), we find 
that the unfolded fields obey the correct chiral commutation relations 

[^a(x),^a’(x')] = inadaa’Sgnix - x'), (23) 

on the whole interval x, x' 6 [0,2L]. In terms of these unfolded fields, our Hamiltonian for the relevant perturbations arising 
from umklapp scattering and superconductivity reads 

J rlL 

dx [gu(-*) cos(2[^r(x) - ^l(x)]) + gsc(-^) cos(^r(x) -h ^z,(x))] 

0 

(24) 


where the position-dependent couplings gu(-*) and gsc(-^) have support on x 6 [0, L\ and x e [L, 2L\ respectively. The unfolded 
system consists of two adjacent regions. Between x - 0 and x - L, we have a region of superconductor where the field 
0(x) = -(^R + ^l)I2 is pinned by the term cos[20] to the value 6* = (n + I / 2 ) 7 r for integer n. From x = L to x = 2L, there is a 
region of Mott insulator where the field (p - (^r— ^l)I2 is pinned by the term cos [40] to the value 0* = 1 /2(m + 1 / 2 ) 7 r for integer 
m. The spectrum is completely gapped, except possibly at the boundaries between the two regions, x = 0 and x - L, where the 
parafermion states we describe emerge. The unfolded system is identical to the topological insulator edge state system studied 
in Ref. Il29l . which tells us that our original system contains Z 4 parafermion state at its ends. 

Let us reproduce the essential parts of the derivation here. We define the total charge and total spin operators for the system, 
according to 


nS = 6(2L) - 0(0), (25) 

nQ =0(2L)-0(O). (26) 

Despite the fact that the fermionic fields must be continuous, the bosonic fields 0 and 0 may jump by integer multiples of 2n, so 
that S and Q can be nonzero in spite of the periodic boundary conditions. The spin S of the system takes integer values, and is 
conserved mod(4) so i = [0,1,2,3) (measured in units of /i/2), whereas the charge Q takes half-integer values and is conserved 
mod(2) so q - (0, j, 1, |) (measured in units of e). Since we may only add integer amounts of electronic charge to our junction, 
we must restrict our value of the charge to be g' e (0,1) ll49ll . The state of our system is then defined by |i, g'). Since every 
physical electron carries one unit of spin, this means that for the charge state q - 0 , only the two total spin states s e ( 0 , 2 ) are 
permissible. Similarly, q - 1 requires s e (1,3). Hence, we have a total fourfold degeneracy of the ground state. 

To see explicitly the parafermionic statistics of the bound states at x = 0 and x - L,we write the pinned values of the fields in 
terms of integer-spectrum operators m, ni, na as 


fll.2 



n 

= 7r«i,2 + 2' 


(27) 
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FIG. 3: Experimental setup for the measurement of the Stt periodic Josephson effect. Two superconductors underneath a Rashba wire are held 
at a phase difference Fractionally charged bound states are indicated by stars. 


These operators then have commutation relations 


[m, n\\ — —, (28) 

n 

[m, 112] - 0. (29) 


The total spin operator is given by 


and the total spin in the system is (5) = (0i - 62)!^ - n2- n\. Then the parafermion states obeying 

T1T2 = 


are given by 




where Tg is the raising operator for charge Tq\s, q) = |i, ^ + 1). 


(30) 


(31) 


(32) 


IV. JOSEPHSON EFFECT & SHAPIRO STEPS. 

A zero-bias conductance peak is a possible experimental signature of localized Majorana fermions. However, such peaks 
could arise from other mechanisms, e.g., disorder ll^l40ll . and do not directly indicate the degeneracy of the states involved. In 
particular, the fourfold degenerate bound state we propose would lead to the same zero-bias anomaly in transport measurements, 
albeit at vanishing magnetic field. To uniquely discriminate these particular bound states, we instead propose to discover their 
presence via the periodicity of the Josephson effect, similar to the corresponding proposal for Majorana fermions IIT0ll2^ . 

To investigate the effect of the zero-energy bound states on the Josephson effect, we follow the logic of Ref ll26l l29l and 
consider an arrangement with two superconducting contacts with phase difference 0sc placed under a Rashba wire partially 
gapped by spin-umklapp scattering (see Fig. |^, in an analogous arrangement to the experimental setup of Ref. Il2^ . The wire 
adjacent to the edges of the superconductors will host zero-energy modes with charge e/2, which will dominate the transport 
at low energies and for a short junction. Tunnelling of a single quasiparticle through the junction changes the parity of the end 
states. In order to satisfy the boundary conditions due to the applied superconducting phase,/our e/2 quasiparticles must tunnel 
via the bound states, leading to an 8;7r periodic Josephson effect ll26l . The time-reversal breaking of the superconducting phase 
difference causes a slight lifting of the fourfold degeneracy, but for realistic parameters this shift is negligible ll26l . 

Several works BTHtSII have suggested that in a spinless, one-dimensional system, the greatest achievable topological degener¬ 
acy is twofold, leading to the statement that only Majorana fermions can exist in one-dimensional systems. Our system does not 
contradict this theorem because in our case the ground state degeneracy is not entirely topological. Indeed, it can be viewed as a 
twofold topological degeneracy combined with a twofold degeneracy due to time-reversal symmetry m. This second degener¬ 
acy can be lifted by local TR symmetry-breaking perturbations such as a magnetic field. In that case, only the topological part of 
the ground state degeneracy survives, and one recovers the 4;7r periodicity of the Josephson effect seen for Majorana bound states 
fTOl l(50l . In real material samples, time reversal symmetry may also be weakly broken by magnetic impurities, thereby lifting 
the fourfold degeneracy to a twofold one, even for a very long wire. In either case, the undriven Josephson current is no longer 
87 r periodic. This raises the question of whether remnants of the 87r periodicity can be observed in such a nonideal setting. 

A possible answer was proposed for similar problems in Majorana nanowires Il44l . where a An periodicity is reduced by 
parity-flipping perturbations to a trivial 2n periodicity: by driving the current in the junction at a finite frequency, we allow 
Landau-Zener tunneling between the different low-lying states. Then, Shapiro step measurements m can still distinguish 
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FIG. 4: Numerical solutions to Eq.(|36) with a d.c. bias current, showing different periodicity behaviour depending on the amplitude of the d.c. 
current. We have made the generic choice of 4 = 10i,„ = lOOip = 1 for the relative size of the contributions to the Josephson current from 
Cooper pairs, Majorana fermions and parafermions. Panel (a) shows the superconducting phase changing in In steps for the choice ao = 1-5, 
whereas panel (b) is plotted for ag = 1.1, and shows residual small steps at In and dominant steps at 4;r which corresponds to the dominant 
transport being via tunneling of Majorana fermions. In panel (c), we choose ag = 1.08, and recover jumps of 8;r in the superconducting phase 
due to tunneling of parafermions. 


higher periodic components even when those are very weak, as in the case of a 47r periodicity recently reported in experiments 
on Majorana bound states Il23ll25l . 

In our finite-length Rashba wire, there will genetically be a nonzero overlap of the modes at each end of the wire and so the 
degeneracy between the modes will be split, although this effect is exponentially suppressed in the length of the wire. However, 
in the case of a driven junction, Landau-Zener tunneling gives us access to all the low energy modes, even when they are subject 
to a small splitting. Allowing for the possibility of Josephson tunneling of Cooper pairs, and for the tunneling of Majorana 
fermions and Z 4 parafermions through the weak link in our system, the total Josephson current flowing is given by 

I((p(t)) = ic sin((^(f)) + im sin(^(f)/2) -t- ip sin(^(f)/4). (33) 

The amplitude ig accounts for the current due to the tunneling of Cooper pairs (this is the critical current above which the junction 
becomes “normal”). The parameters i„, and ip similarly account for the tunneling of Majorana fermions and Z 4 parafermions 
through the junction. The Josephson equation relating the rate of change of the superconducting phase tp to the voltage across 
the junction reads 

m = ^v(f). (34) 

n 

We current-bias our Josephson junction using a constant current with a small a.c. component Iq -h Ii sin(wf). The current 
through the junction consists of two parallel components, a tunnel current given by I{tp{t)) and a resistive current due to ohmic 
quasiparticle transport in the junction. Equating the sum of these two contributions to the bias current, the gauge-invariant phase 
(p{t) is described by the dynamical equation 

n 

lo + h sin(mf) = ic sin(ip(t)) + i,„ sin(ip(t)l2) + ip sin{(p(t)/4-) -H :^ip(t). (35) 

2eR 

Note that these contributions to the Josephson current are periodic under shifts (p(t) (p{t) ■+■ 2n, (p(t) p{t) -H 47r and (p{t) 

ip{t) ■¥ 2>n respectively. Writing the equation ( |35| l in terms of rescaled variables t = 2eRictlfi, and d> = tia>/2eRic, we find the 
equation 


<P(t) - aQ + a\ sin(a;T) - sin(^(T)) - Um sin(^(T)/2) - Qp sin(^(T)/4) (36) 

This equation must be solved numerically. 

In order to see that the Stt periodicity can win out, even when = 1 > o'™ > ctp, we first solve ( [36] l without the a.c. drive 
current, i.e. cri = 0 (see Fig. |^. For a generic choice of the d.c. bias current, ao, the superconducting phase is 27r periodic, 
reflecting the dominance of the Cooper pair tunneling (Fig. |^). However, we find that as we approach a critical value of ao, 
we first see the An periodic term due to the Majorana contribution (Fig. [^) and then the Stt contribution from the parafermions 
(Fig.|^) become dominant. 

The winding up of the superconducting phase is not a property which is easy to directly measure, so to give an experimentally 
accessible measurement, we must drive the Landau-Zener transitions which will allow us to see the degeneracy of the low-lying 
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FIG. 5: A plot of the current-voltage relationship for the system driven with a small amplitude a.c. current with a large d.c. offset (corresponding 
to Offl = 1.08 -I- Sag and ori = 0.01). The excess current Sag over the critical value at w = 0 is plotted on the vertical axis. The Josephson 
junction has a dominant Itt periodic component, and sub-leading An and 8 ;r periodic components (4 = lOi^ = IOO 4 = 1). The dashed curve 
shows the system driven at a frequency (jj\ = 0.04 (comparable to that of the 2n periodic component of the Josephson current), and shows 
all integer Shapiro steps. The solid line is the system driven with an a.c. component with frequency (1)2 = 0.01, and only even Shapiro steps 
survive. Finally, the dashed-dotted line shows the Cji = 0.00052 behaviour, and the disappearance of all the Shapiro steps which are not a 
multiple of 4. Note that on the same scale as the Wi and 0)2 plots, the W 3 shows very small variation, so its amplitude has been increased by a 
factor of 7 to clearly show the disappearance of the Shapiro steps. 


modes. To do this, we switch on the small a.c. component to the bias current ai (experimentally achieved by irradtiating the 
junction with microwaves). Tuning the driving frequency allows us to access three distinct regimes, in which steps occur in the 
experimentally-measurable I-V curves for the junction at distinct multiples of the driving frequency d). For high frequencies, 
we hnd Shapiro steps are present at all integer multiples of the driving frequency, indicating that the transport is dominated by 
conventional Josephson tunneling of Cooper pairs through the junction. Reducing the frequency of the drive, we recover the 
results of Ref. El, that only the o^f£/-numbered Shapiro steps remain. Finally, for very low frequencies, there exists a regime 
in which only every (4n + 1)'^ Shapiro step survives (see Fig. |^. Whilst the appearance of extra Shapiro steps can be caused 
e.g. by disorder, alternative mechanisms to the one described by which Shapiro steps may disappear seem to be unknown. As 
the disappearing Shapiro steps are robust even when the 2n and An contributions to the Josephson current are dominant over the 
Stt component, they therefore provide a highly selective test of the existence of Z 4 parafermions. 


V. DISCUSSION. 


To summarize, we predict the existence of fourfold degenerate zero-energy bound states protected by time-reversal symmetry 
when inducing superconductivity in strongly interacting Rashba wires. The bound state operators acting in the ground state 
manifold satisfy a Z 4 parafermionic algebra, and transitions between the degenerate ground states occur via the tunneling of 
fractional charges e/2. We propose an experimental scheme allowing us to observe an Stt periodic Josephson current and 
associated Shapiro step structure, which would be a dehnitive signature of our predicted bound states. 

The opening of the spin-umklapp gap at 4 = 0 is rather generic. In a wire with RSOC and electron-electron interactions this 
gap will lead to a reduction of the normal-state conductance of from 2e^ Ih to jh as the chemical potential approaches the Dirac 
point. A similar phenomenon has been observed in GaAs wires ll46ll and our results may provide an alternative interpretation for 
this experiment if RSOC is sufficiently strong in the semiconductor nanowires used. A conductance measurement as a function 
of the chemical potential in InSb or InAs wires, where RSOC is typically stronger, could demonstrate the proposed reduction of 
the conductance conclusively. We would like to point out that the required strong interactions have already been seen in these 
wires El. Umklapp scattering has also recently been invoked to explain the observed conductance reduction in InAs/GaSb 
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topological insulator edge states ll48l . Introducing a weak superconducting proximity effect in either of these types of wires 
or edge states will then lead to the creation of bound states and allow the observation of the Stt periodic Josephson current 
component. 
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Appendix A: Interactions 

The interaction Hamiltonian for a density-density interaction reads 

^ ^ Xj f -r2|)iA;r2(r2)iA<r,(ri). (Al) 

0-10-2 


and the operator lAl-(r) creates a particle with spin cr e {"f, i) at position r = (jc.y). We transform this term using the Schrieffer- 
Wolff transformation, and then project onto the lowest subband. Due to translation invariance, the momentum in x direction 
remains a good quantum number. The Hamiltonian in momentum space in this sector is 

H^Ho + Vp + Vsf + V.X. (A2) 


where the kinetic energy takes the form 


/To = 


a-,p 


— -I- aRcrp 

Zm 


^p,a-^ P,a-^ 


(A3) 


and the operators and 4>p,cr create and annihilate an electron of spin cr and momentum p. The potential terms are then Vp, 
and the density-density interaction becomes 

a-10-2 p,p',q 


where we have introduced V{q,y) as the partial Fourier transform in x direction of the physical interaction potential y(|r|) s 
V(x,y), 


V(q,y) - J dxe '^^y(x,y) 

and defined V{q) = V(q,y - 0). Next, the spin flip and spin exchange terms are given by 

Vsf = i Zo- l 2 p,p',q U(q)(2p' - q)(2p + q)fpl+q,a-^l,_q^^^p'-a->Pp,-cT, 

Kx = Z U(q){2p' - q){2p + 

with an effective potential U{q) given by 

mat -Air” 2 2 fel-Z27 2/7 ~ / Zr \ 

U{q) = 2,12 3 / , , - I dzidzidz^e ^ //„(zi - Z2)//i(zi)//i(z2)e Hn{zi)V\q,——\, 


(A5) 


(A6) 

(A7) 


(A8) 


where N„(z) denote Hermite polynomials which emerge because they are the transversal eigenfunctions due to the harmonic 
confinement in y direction. 

We now project these terms onto the linearized low-energy states of the system. For the density-density term, we obtain 

yp = y( 0 )^ I dxpa,Ax)Pa,Jx) 

a,a- ^ 


(A9) 
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with the fermionic densities Paa-ix) - if/l.a-(x)i//a(r(,x). Due to momentum conservation in x direction, the spin flip terms may only 
mix terms near q = 0, so Vsf reads 


Vsf = - 2 U( 0 )J' +h.c.) (AlO) 

Finally, the spin exchange terms only lead to a small change in terms already present in Vp, so we can account for them by a 
change of the interaction parameters, but we should keep separate the one term which may not be written as a density-density 
interaction; 


Vs = 2k\U(kF) J" 

We now bosonize these terms. Using the standard identity that = dx4>a,i7l2n, we find that at weak coupling, the complete 
interaction Hamiltonian V — Vp + Vsf + Vsx becomes 

^ =0/ + + ^ f dx^l-8nO)-2k^pU(kF)jd,,pA<P-+2k^^U(kF)d.M^ 

f f dx (e‘(2«+-20-) + jj.c.) (All) 


We note at this point that the Luttinger parameters for the two species are equal in the weak coupling limit. The coupling 
between the species given by the terms K^dx^+djc(f>- where - -8U(0) - 2k^pU(kf) and Kod^O+djcd- where Kg = 2kj^U(kf). 
These can be removed by going to the charge-spin basis described in the text above Eqn.(lO), which yields Eqn.(lO) with new 
Luttinger parameters 


J^2 _ K+Kgjl 

~ 1IK+k^I2^ 
j ^2 _ K-K0I2 

^cr \IK-K4,I2^ 


(A12) 


and renormalised Eermi velocities Vp - VfKp and Vq- = VfKg-. Vf - ctr is the Eermi velocity for the non-interacting modes. 
Note that these expressions are only valid at weak coupling. However, the division into a pair of non-interacting Luttinger 
Hamiltonians with three competing cosine interaction terms forms our prototypical model for the strongly interacting case we 
consider. 


Appendix B: Renormalization Group 

We calculate the scaling dimensions of the operators above, and also the second order RG flow for the parameters Kp and Kcr- 
We diagonalize the non-interacting Hamiltonian (Eq.(lO) in the main text) by introducing the fields ipaa - a yfK^j(pa - 0al 
where a - p,cr, and as before a - R,L, which gives 


Ho^ 


2jt 


J" dx [{dxPpRf + {dx^Pplf) + ^ {{dxV’a-Rf + (dxV’a-Lf) ■ 


(Bl) 


The scaling dimensions of the operators are calculated by normal-ordering them with respect to the creation and annihilation 
operators of the diagonal Hamiltonian (Bl 1 . We find 


Vs = 


gs / L 


{2na)^ \2nal 


f [2 yfioA : ■ 


(B2) 


Asserting that this term cannot change as a result of the RG step a ^ a —da - a(\ + d{) and gs gs + dgs gives the first-order 
RG equation for gs as 


dgs 

dt 


2 -A 


gs- 


(B3) 


Under normal ordering and point splitting, the umklapp term becomes 


Vu = 


(27ia)^ \2na) 


-2(K„+K^) 


J" dx : cos [2 V2(0p - 0a-)] :, 


(B4) 
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and so has the first-order RG equation 


dgv 

d€ 




The superconducting term has scaling behaviour 

gsc 


Tsc - 


(2na)^ 


" dx[: cos[ V2(6»p H- 0o-] : H- : cos[ V2(0p - 0a-] : j. 


giving an RG equation 


dgsc 

M 


2 - 


2Ka 2K„ 


gsc- 


Continuing to second order in the spin density wave term, we find the real time partition function takes the form 


= ) 
8(27ra)"^ V 2m I 


L 


/ 


dx\ dx2 dti dt2 : cos[2 V20o-(xi, ti)] :: cos [2 ^0a(X2, t 2 )] : . 


(B5) 


(B6) 


(B7) 


(B8) 


We then normal order the two individually normal ordered cosines, and keep only the part which gives us a new scaling compared 
to the first-order term to get 


Z 2 = - 


8(2mY 


I 


dxi dX2 dt\ dt2 [maMa] : cos[2 ^ 6 a{x\,t{)\ cos [2 V20cr(x2, fa)] 


(B9) 


where we have defined m* - a-iva(t2 -fi)± i{x2 - xi). Because they are now under normal ordering, we can expand the cosines 
for small X 2 - xi, to get 


Z 2 = - 


gj^ 

(27rfl)4 


/ 


dxi dX 2 dti dt 2 [maMa] '^^^"{X2 - X\f : {dx9a) 


^2 . 


(BIO) 


This is a renormalization of the coefficient of {dx9aY, that is to say a renormalization of Ka- We now re-express this as a term in 
the first-order partition function of the operator {dx9a)^\ we therefore shift to centre of mass x - xi + X2 and relative x - X2- x\ 
coordinates (same for t), and leave the integration over the centre of mass coordinates alone. In order to calculate the coefficient 
of the term (dxOaY, we do the RG step a ^ a — da - a{\ + df) inside the integral only, which gives us an integral form for the 
coefficient ki 






- iVatf + 

{2Ti)^Va 


which turns up in the renormalization of Ka, dKa = K\g^d(. We can compute this integral exactly, to get 


(BID 


1 ^P:T[2|Ka - 3/2] 

(27r)4v, 2Y[2IKa] 


(B12) 


This term has a weak, linear dependence on Ka, so only allowing for small changes of Ka we may treat it as approximately 
constant. 

Similar calculations for the terms in gu and gsc give us the coupled RG equations below 


dKg- ^ -y yj -y 

(B13) 

dKp ^ 

— ^ -K2g^Kp + K2g^^. 

(B14) 


/C 2 is given by the integral 




X2 = 


(2mr 


r r“ i 1 

dt dxx^da\[(a - ivpt)^ + x^y^^" [(a - ivat)^ -(- |, 

\ J—CC % J—CO '■ ^ 


(B15) 
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and by 


K3 = 


2a 


(i + sl) 


1+1 


{2nar 


X O poo . ^ . 

dt I dxx^dal[(a - ivpt)^ + [(a - iva-t)^ + 

oo %J —oo '■ ' 


(B16) 


Note that the use of the full, second order RG equations does not change qualitatively the result given in the main text using 
only the first order RG equations, as the second order terms only result in small changes in Kp and Ka-. 
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